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SUMMARY 



Charts are presented for the calculation of the crit- 
ical stress for local instaljility of col-umns with I—, 
channel, and rectangular— tube section. These charts are 
intended to replace the less complete charts published in 
NACA Technical Note No, 743, The values used in extend- 
ing the charts are computed "by moment— distribut ion methods 
that give somewhat more accurate values than the energy 
method previously used and also make it possible to deter^- 
mlne theoretically which element of the cross section is 
primarily responsible for instability. 

An experimental curve is included for use in taking 
into account the effect of stresses above the elastic 
range on the modulus of elasticity of 24S-T aluminum alloy, 

A determination of the dimensions of a thin— metal 
column for maximum critical stress with certain given con- 
ditions is presented. 



INTRODUCTION 



One of the important requirements in the design of 
thin-metal columns for aircraft is the determination of 
the critical compressive stress at which local instability 
occurs. Local instability of a column is defined as any 
type of instability in which the cross sections are dis- 
torted in their own planes but are not translated or ro- 
tated. 



2 



The critical stress for local instability can 
usually be given in terms of the geometry of the sec- 
tion, the properties of the material, and a coefficient, 
Reference 1 presented, charts for the determination of 
such coefficients for columns of I—, channel, and 

rectangular— tube sections, These charts, however, con- 
tained relatively few curves and in some cases required 
interpolation over a wide range. 

In order to make the charts of reference 1 more 
nearly complete and to reduce the necessary range of 
interpolation, each chart has been extended to include 
eight intermediate curves, The values used in extend- 
ing the charts are computed by moment-distribution 
methods that give somewhat more accurate values than the 
energy method previously used and also make it possible 
to determine theoretically which element of the cross 
section Is primarily aresponsibl© for instability. 

The present report includes the extended charts, 
aifing with tables of the values used in preparing the 
charts, and is intended to supersede reference 1. An 
experimental curve is included for use in t/aiing into 
account the effect of stresses above the elastic range 
on the modulus of elasticity of 24S— T aluminum alloy« 
A determination of the dimensions of a thin— metal colxuan 
for maximum critical stress with certain given conditions 
is presented. 



SYMBOLS 



A 



cross— sect ional area 



width of end or narrower wall of rectangular tube or 
of plate element of I—, Z— , or channel section 



effective flexural stiffness of plate per unit length 
i 3 _ 

r T^Et 1 




modulus of elasticity 



h 



width of side or wider wall of rectangular tube 
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k nondimensional coefficient dependent upon relative 
dimensions of cross section 

kgQQ section coefficient 

t thickness 
III 

S stiffness in moment— dl stri"bution analysis for far 
edge free (no support and no restraint against 
rotat ion) 

S^^ stiffness in momentr-distr ibuti on analysis, for £ar 
edge supported and subjected to sinusoidally 
distributed moment equal and opposite to moment 

applied at near edge 

e restraint coefficient, a measure of relative resis- 
tance to rotation of restraining element at edge 
of plate 

\ half wave length of buckle 

\x Foisson's ratio 

Oqj. critical compressive stress 

n nondimensional coefficient that takes into account 
reduction of modulus of elasticity for stresses 
above the elastic range. Within the elastic 
range, n = 1. 

Subscripts: 

I* flange 

V web 

b end or narrower wall of rectangular tube 
h side -or wider wall of rectangular tube 

POHMUIAS yOE CEITICAL STHESS 



Por an I—, Z— , or channel soction, either of two 
formulas given in reference 1 , may be used for calculating 
the critical compressive stress. The two formulas are 



and 



'cr 
r\ 



12(1 



(2) 



The corresponding formula for a rectangular-tuTae section 
is given in reference 1 as 

f cr _ ^^^^^h' 
i\ ~ 12(1 - |j.«)h» 

In using formulas (l), (2), and (3) when the stresses 
are ahove the elastic range, CTcr/ti is first evaluated, 
and Oqj, is determined from this value hy means of the 
curve of figure 6, The relationship between and 
^otI"^ will he further discussed in another section of 
this report, 

DlSCUSSIOff or CHARTS 



All of the quantities on the right-hand side. of 
equation (l), (2), or (s) are knovrn except the value of 
the coefficient k^, kp, or k. This value may be read 
from the appropriate chart (figs, 1 to 5) after the nec~ 
essary dimension ratios are computed and applies v/henever 
the length of the column is greater than several (3 or 4) 
times the width of the widest plate element^''. 

In general, when a column of I-, channel, or 

rectangular— tube, section fails by local instability, one 
of the two elements (tveb and flange or end wall and side 
wall) of the cross section may be said to be primarily 
responsible for the instability; that is, as the load 
approaches its critical value, this one element is no 



5 



longer capalile in itself of supporting the loads imposed 
oh it ir/ithout "buckling and ireitiuires a ce'irta of 
restraint from the other element of the cross section in 
order to delay tuckling until that load for which the 
cross section as a whole 1)600068 unstable is reached. 
The charts show which element of the cross section is "be- 
ing restrained against "buckling hy the other element, A 
dashed line is drawn on each of the charts (figs. 1 to 5) 
connecting the points for which the two elements are 
equally responsible for the instability of the section, 
This line divides the chart into two regions: In one 
region the web (or side wall) is primarily responsible 
for instability and in the other region the flange (or 
end wall) is priiaarily responsible for instability, A 
column with a given cross section will fall into one of 
these two regions, dei^endiag on the values of the various 
dimension ratios. 



RELAOJIONSHIP BETWIEIT AHD acr/n 



Figure 6 shows the relat i onsfcip • between a^j. and 
CTgp/Ti as determined from tests of 24S-T aluminum— alloy 

columns of Z— , H— , and channel section, either formed from 
flat sheet or extruded. This figure was prepared by plot- 
ting the experimentally determined values of a^j. as 
ordinates against the values of cr^^/n as abscissas. The 
values of agp/ri were computed according to equation (1) 

and the chart of f igure - 1. o,r- 3. The results of the tests 
are dirscussed in more detail in reference' 2, 

Similar experimental data for materials other than 
84&-T aluminum alloy are not now available, and further 
study of this subject seems desirable, ■ 



METHOD OF PHEPARIiTG CHARTS 



Values of k^, ky, and k used in preparing the 
charts (figs. 1 to 5) were computed by an application of 
the principles of moment distribution to the stability of 
thin plates. This method is presented in detail and one 
example of its application is given in reference 3. 



Aq alternate procedure, which makes use of the 
charts presented in references h and 5» was used in com- 
puting scane of the values for the charts of figures 1 to 
3' An example of the application of this method follows: 

It is desired to determine the value of the coeffi- 
cient fciif for a column of the cross-sectional dimensions 
shovm in figure J. It is necessary to predict which ele- 
ment of the cross section will be primarily responsible 
for the ir'stability. The calculations will then show 
whether this prediction is correct. For a section with 
a relatively wide flange, such as that shown in figure 7, 
the flange will probably be primarily responsible for 
failure. On this assumption, the detailed procedure for 
determining ki^^ or kp is as follows*. 

1, Compute the ratios trj/t^Tj and byy/bp. 

2* Assume a value of V'^'j where X is the half 
wave length of the buckle. 

3« Compute A./bp from tha equation '^JT 

Assume several values of ky^j-. In order to avoid 
the necessity for interpolation in the tables of reference 
6, the values of k should be assumed for the part that 
is not primarily resoonsible for instability because these 
values are the ones that must be used to enter the tables. 
Unless previous exnerience has revealed the approximate 
r^nge of such values to be assumed, this range may be 
estimated by the use of the relation between the values 
of k for the two parts of the cross section given in step 5 
and by the fact that the value of K for the part which is 
primarily responsible for instability will be somewhere 
between that for simple support and that for fixed edges. 

5. For each value of "k^p compute kp from the 
equation 



which was obtained from equations (l) and (2) and the 
..assumptions hat sttefs ig .uniform across ihes.ect ion, 

6. Using the assumed values of \/hy and k^, 
evaluate the quantity s'\/il3/^)^ from the tables- of 
reference 6, where 







12(1 - p,*) 

o 4S^\bs, 
7, Compute € = (gee. reference 4), where 



The formula is 



riiltj,=* 



12(1 - ti^) 



^t'r. the values of e from step 7, determine 

from the chart of figure 3, reference 4! 

9. Plot kj, from step 5 and kj from step 8 as 

?hp'?^i^ against either of the two, values as abscissa. ' 

value nl T^r ^iJ""^ ^^^^^ °°"«°t 

value of kp for the particular value of \./bj. 

of jt/by Sepeat steps 2 to 9, assuming different values 

* . • V. • 1 .... ' 

x/-h ^^%v,^^°- values of kj, from step 9 against 

The minimum of this cur^e, gives the required value 
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iE- V 

If the calculations indicate that S v is nega- 
tive, the prediction that the flange is the primary 
cause of instability is wrong. In such a case, the 

calculation must "be carried out with S^^^y instead of 
S^^Y^ in step 6, and with the chart of figure 3, refer- 
ence 5, in step 8, In addition, all the subscripts I 
v/ill "become W, and vice versa, 

Ihe results of the procedure outlined herein as 
applied to the problem of figure 7 are given in table I, 
The values of k^i in the last column of table I were 
determined according to step 9, If these values of kji 
are plotted against X/by, the minimum value is found 
to be about 0.73. The value of can be computed from 

the formula given in step 5. 

niables II to VI give the minimum values of k^, kp, 

and k used in the preparation of figures 1 to 5, All 
of the values of k and k^^ in these tables except those 
marked a were computed either by the method just out- 
lined or by the moment— distribution method discussed in 
reference 3. The values of ky were then computed by 
the eqLuation given in step 5, The values marked a are 
those computed by the energy method and used in the prep- 
aration of the charts of reference 1, 



DIMENSIOirS OP THIIT-METAL COLUkKS POE 
MAXIMUM CEITICAL STHESS 



Eq.uation (l) gives the critical stress for an I—, 
Z— , or channel column in terms of the width and the thick- 
ness of the v;eb. The effect of the presence of flanges is 
taken into account in the evaluation of the coefficient 
k^, For the purpose of studying the dimensions that give 
maximum critical stress, the foriri of equation (l) is pre- 
served but the concept of certain terms is generalized. 

The ratio b/t of a plate may be called the aspect 
ratio of the plate, A corresponding quantity that ex- 
presses the "section aspect ratio" for a thin— metal col- 
umn is the area of the section divided by the square of 
some thickness. If, therefore, equation (l) is written 



9. 



(4) 

n 12(1 - n^) V A y 

then the value of the section coefficient kge^ is a 
measure of the effect of the shape of the section hy/hv 
on Oq^/'^ for a given section aspect ratio A/t^^ and 
a given value of t^/tp. In order to show that kgec ^® 
dependent on only hj,/ by^ and ty^/tp, equation (1) is 
set equal to equation (4), with the result that 



k 



From the geometry of the section ( Z— or channel), 

A = h^/tyf + Shptg. If this value of A is substituted 

in equation (5) and the equation is solved for kgeci 
the result is.: - 



*-sec 



The value of depends on only • hp/b^ and t^/ty, ■ 

and the value of kgec therefore also depends on only 
these two ratios. 

In figure 8 the values of ^^qq as determined by 
equation (6) are plotted for channel-^ and Z-section" col- 
umns, and in figure 9 similar values are plotted for I- 
section columns, 

A method exactly analogous to the foregoing method 
can be applied to rectangular tubes. In this case, equa- 
t ian (4) is writ t en . . : : 

. gcr _ ^ ' Tf^E /*h"^>,° ■ - . 
^scc TT- r- ( j (7) 
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and the formula for kgec 'becojnes 

la figure 10 are plotted the values of kgeo for rec- 
tangular tuTjes, ad determined "by equation (8), 

As a practical problem in the determination of the 
dimensions of a thin— metal column for the development of 
maximum critical stress, consider a flat strip of metal 
of constant thickness which is to "be formed into a ^ or 
channel section. In the formed section, tw/^F ~ 1« '^'^^ 
section aspect ratio A/tw^ is equal to the v/idth of this 

strip, or the developed length of the final cross section, 
divided hy the thickness,. When bent to form a channel— or 
Z-section column, this strip of metal of constant thick- 
ness develops the highest or^p/'n for local instahility if 
the hends are so located that the ratio of flange width to 
web width 'bj/t'^ is equal to about 0,41, which is the 
maximum of the curve for ^y/tji =1 in figure 8, 

Regardless of the thickness used in the definition 
of the section aspect ratio, the maximum -ivalue of Cfcr/T) 
for a given value of the section aspect ratio will occur 
at the same values of by/b'fif for a particular value of 
*w/*P» maximum for each t'jf/tp ratio therefore re- 

veals the shape — that is, the value of bj/b^ — that the 
I—, Z—, or channel section should have if maximum o'cr/'H 

is desired. The same reasoning holds for the rectangular 
tube, (See fig, 10.) 

Equations (l) to (s) and figures 1 to 5. are probably 
more useful to practical designers than the more general 
equations (4) and (?) and figures 8 to 10^ The curves of 
figures 1, 3, and 5 have' theref or e been redrawn in figures 
11 to 13 with dashed lines added to show the percentage of 
the maximum value of o-Qj./n that can be developed for given 
values of t^/tj. and A/t-j/^ when bp/b^ is varied. The 

position of these lines is independent of the thickness 
used in the definition of the section aspect ratio. It is 
of interest to observe,, by comparison of figures 2, 4, and 
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5 with figures 11 to 13, that the line of maximum values 
tears no apparent relation to the line that shows the 
dimension ratios for which the web and flange (or end 
wall and side wall) are equally responsiljle for the in~ 
stal)ility of the section. 



C OUCIUSI OLIS 



1, The critical compressive stress at which cross- 
sec.ional distortion begins in a thin-wall column of I-. 
f* or channel section is given hy either of the follow- 
ing formulas: 



12(1 _ |a^-)V 



or 



^cr _ kpTT^Etj,* 



12(1 - ti^)hj^ 



where 



"by width of weh 

hj. half width of flange 
flange for Z- and 

and kp nondimensi onal 
appropriate 



for I-section, total width of 
channel section 

coefficients read from the 
char t 



E and p. Youngis modulus and Poissonis ratio for the 
material, respectively 

t^ and tj thickness of web and flange, respectively 

Tf» nondimensional coefficient that takes into 

acccuat reduction of modulus of elasticity 
lor stresses above the elastic range. Ifithin 
the elastic range, n = 1, 
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For a rectangular— tuTie section 



^cr , ^ ^^h^ 
^ 13(1 - |j,^)h^ 



where 

k nondimensional coefficient read from appropriate 
chart 

h and width and thickness, respectively, of side or 

wider wall of rectangular tube 

2, For stresses above the elastic range, the critical 
compressive stress is determined from a curve that gives 
the relationship between a^^ and Oqj,/t\ for 24S— T alur- 
minum alloy. 

3, The charts of values of k are divided into two 
regions: In one region the web or side wall is primarily 
responsible for instability and iu the other region the 
flange or end wall is primarily responsible for instability. 



4, The equations for critical stress are also pres— " 
ented in general form with the ratio b/t replaced by the 
section aspect ratio A/t^, where b is the width and t 
the thickness of an element of the cross section, and A 
is the area of the cross section. Prom these general eqtua— 
tions, charts have been prepared that reveal the effect of 
shape alone on the critical stress for local instability. 
The shapes that give maximum critical stress bear no appar- 
ent relation to the proportions for x\rhich the web and flange 
(or end wall and side wall) are equally responsible for the 
instability of the section. 



Langley Memorial Aeronautical laboratory, 

National Advisory Ccm.'riittee for Aeronautics, 
Langley SMeld, 7a, 
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TABLE II 

CALCULATED MINIMUM VALUES OF ky^ FOR I-SBCTIONS 





0.5 


0.6 


0.7 


0.8 


0.9 


1.0 


1.2 


1.1+ 


1.6 


1.8 


2.0 




^6.05 
^6.1^6 
^6.61 
^6.68 
*6.72 
6.72 

"67711 










^J+.oo 

*il..82 

as. 01 
1^.86 
i4..62 
U.12 

5.58 










nl+'OO 

^1;.06 

"l+.og 
^5.89 

5.57 
2.97 
2.25 
1.72 


0 

.05 

.10 

.15 

.20 
.25 
.50 

.i^75 

.50 

.525 


































































646 

6.1^9 
6.51 
6.52 
5.46 
5.96 
5-51 


6.12 
6.15 

6.09 

5.7t 
5.01). 


5.71 
5.59 

B 


5.26 

1^.07 
5.52 


U.32 

5.1+9 
2.97 


5.99 
5.59 
5.07 
2.55 


5.7? 

5-54 
2.78 

2.25 


3.66 

5.17 
2.5L 

1.98 
















6.75 
6.75 
6.59 
6.21 

5.95 
5.54 
11.25 
5.50 
2.19 


1^.51 




5.02 


2.61 


2.12 


1.76 


1.50 


1.50 


1.10 






















.55 






















.56 

.575 

.60 

!^0 
1.00 






















l,.if 
5.12 

2.1+7 
1.65 


5.22 

2.1^6 

I.9U 
1.29 


2.65 
2.02 
1.58 
1.07 


2.25 
1.71 
1.55 
.90 


1,18 

.78 


1.51+ 
1.17 

•f^ 
.61 


1.26 
.96 

.1+9 


.ho 


.90 
.66 

.51 
.55 


.76 

"1 



■computed by energy solution (reference 1). 
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TABLE III 

CALCULATED MINIMUM VALUES OP kp PGR I-SECTIOKS 
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0.5 


0.6 


0.7 


0.8 


0.9 


1.0 


1.2 


1.4 


1.6 


1.8 


2.C 


0.0 
.2 


ai.288 










ai.288 
ai.019 










ai.238 


















a. 623 


















ai.193 












a. 852 










.8 

1.0 
1.25 

l.lj.29 

1.667 

1.786 
1.818 

1.905 
2.000 
2.105 
2.222 
2.500 
2.857 
5.333 
I|..000 

5.000 
6.667 


a.567 
.547 
.527 


















ai.l47 
1.120 
1.119 


.550 
.532 


0.634 
.607 
.590 
.568 


0.685 
.646 
.633 
.605 


0.729 

.678 
.651 


.777 
.752 
.725 
.703 


0.883 
.840 
.825 
.798 


0.960 

.937 
.920 
.886 


1.029 

i.oo4 

.982 
.985 


1.07S 
1.064 

i!o^ 


.518 

•^^^ 

.1^87 
.1,84 


1.096 






























































.lu95 
.f< 
.471 
.575 

.287 
.210 

.145 


.528 


.567 


.611 


.655 


.762 


.862 


.962 


1.049 


1.098 




.500 
.450 
.566 
.271 
.187 


.485 
.416 
.322 
.228 


.577 
.520 
,460 
.568 
.266 














.269 
.151 

a. 067 
a. 037 


.573 

■505 

.4ifa 
.304 
' .200 

a. 112 


.685 
.615 
.511 
.380 


.799 
.737 
.633 
.409 


.920 
.871 
.J69 
.606 


1 ,025 
1.007 
.925 
.741 


1.101 
1.095 
1.070 
.892 

a. 623 
a.Jcu- 



































acomputed by energy solution (reference 1). 

TABLE iv" 

CALCUIiATED MINIMUM VALUES OP % FOR CHyvNWEL iiND Z-SECTI0i;S 



NtwAp 


0.5 


0.6 


0.7 


0.0 


0.9 


1.0 


1.2 


1.4 


1.6 


1.8 


2.0 


0 

.050 
.100 

.130 
.167 

.179 
.192 

.200 
.227 
.250 
.500 
.350 
.400 
.450 

.475 
.500 

.525 
.550 
.560 

.575 
.600 
.700 
.000 

1.000 


al+.OO 
as. 1^6 
a6.02 
a6.19 
a6.31 










a^.OO 
ail.. 26 
^4.45 










al(..00 
a4.03 

%.o4 


























































^4.50 

^4. So 

4.59 

t 8 
4i^i 
4.11 
3.74 










^4.00 
a5.98 
|5.97 
3.92 

3.76. 
3.26 

i'.^7 


















a6.50 


































^6.43 
6.50 

6.53 


















6.05 
6.10 
6.13 
6.10 
6.06 
5. 06 
5.46 


5.59 
5.59 
5.55 
5.36 
4.95 


5.15 
p. 12 
4.92 
4.62 

i|.l9 


4.79 
4.69 

4.46 
4.09 
5.70 


4.19 
5.98 
3.65 
3.25 


4.00 
3.76 
3.56 
2.89 


3.88 
3.59 
3.15 
2.56 


3.81 
5.48 
2.78 
2.18 
















6.54 
6!^6 

6.4o 
6.07 


4.45 


3.74 


5.26 


2.90 


2.4o 


2.02 


1.69 


1.42 


1.19 






























































S.70 
4.45 
3.54 

2.59 


4.31 

2.60 
l.Ol 


5-46 

'■> ■ ' 


2.91 
2.27 

1.21 


2.49 
1.93 
1.53 
1.03 


2.20 
1.70 
1.56 
.89 


1.76 
1.32 
1.05 

.66 


i.iOt 
1.09 
.04 

.55 


1.19 

.104 


.98 
.72 
.55 
.36 

*-. 


.85 

""45 
.50 



^Cc-i'i^uted. by energy solutloo (refsranoe 1). 
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TABLE V 

CALCULATED MINIMUl* VALUES OP- POH CHANNEL AND Z-SECTI0N3 





0.5 


0.6 


0.7 


0.8 


0.9 


1.0 


1.2 


1.1; 


1.6 


1.8 


2.0 


0 


^1.288 










l.tiOO 
Sil -111 












.2 


















a. 695 


















1.23i) 


■} 










a. 962 










.8 


a. 621 
.598 

•III 
.545 

.513 
.502 
.501 
.iJ-89 

.14.09 


















ai.20lj 
1.19< 
1.15< 


1.0 

1.250 

1429 
1.667 

1.759 
1.786 
1.818 

I. 905 
2.000 
2.105 
2.222 
2.500 
2.857 
3.555 

II. .000 

ij..l|.00- 
I4..800 
5.200 
5.600 
6.000 


0.650 
.617 

.592 
.559 


0.706 

.676 

.655 
.611 


0.772 
.726 
.712 
.671 


0.836 
.791 
• 767 
.727 


.890 
.870 
.852 
.792 


0.98s 
.96k 

.95i^ 
.911 


1.074 
1.056 

1.045 


1.134 
1.110 

x*110 

1.099 


1.161 
1.140 
1.143 
1.157 


1.015 


1.195 






























































• i4-75 

.551 
.270 
.198 
.156 


Cl, 2 

• 5M-5 


.590 


.657 


.725 




.990 




1.11X1 


1.185 




.L89 

.555 
.247 
.171 


■M 

.295 
.206 


.606 














.261 


.598 
.503 

.286 


.515 
.377 


.906 
.807 
.662 

• 489 




1.151 
1.104 
1.013 
.772 


1.1914 
1.18^ 

i.m 
.959 

a. 681 
a. 587 

a.5oe 


,li4.6 


'3U2 


.621 


a. 085 


a. 236 
a. 199 
a. 170 
a. 11^6 
a. 127 


















a. 059 


































a.OUii. 





































aconrputed by energy solution (reference 1). 

TABLE VI 



CALCULATED MINIMUM VALUES OP k POR RECTANGULAR TUBES 



b/h\ 


0.5 


0.6 


0.7 


0.8 


0.9 


1.0 


1.2 


1.4 


1.6 


1.8 


2.0 


0 

.050 
.075 
.100 

.125 

.200 
.300 
.400 
.■500 
.600 
.650 
.661 
.670 
,700 
.750 
.780 

.Soo 

.820 
.850 
.900 
1.000 


a7.01 










a7.01 

^6.45 










a?. 01 


aJ.'aS 




















^4.72 










^6.09 










a6.85 


a4.62 

4.ii4 
4.30 
4.21 
4.12 
3.94 
5.82 

5-7^ 
5.68 

5.38 


















4.68 


4.94 


5.21 


5-42 


5.66 


6.02 


6.51 


6.48 


6.62 


6.71 


4.41 


4.62 


4.82 


5.05 


5-50 


5.72 


6.05 


6.29 


6.45 

6.41 


6.59 

6.57 

6.56 


4.12 


4.30 


4.52 


4.77 


5.03 


5-52 


5.92 


6.20 










































3.86 
5.66 


4.08 














6.40 


6.55 


































2.56 


5.50 
5.42 
3.27 
3.05 
2.75 
2.23 


5.75 


4.04 


4.56 


4.64 


5.28 


5.79 


6.15 


6.59 


6.5I; 
























5.29 
2.75 


5.69 

5.22 












6.39 
6.59 


'a 


1.64 


3.61 


4.00 


4.81 


5.56 


6.10 



Computed by energy solution (reference 1). 
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